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In a medium that contains a neutrino background in addition to the matter particles, the neutrinos 
contribute to the photon self-energy as a result of the effective electromagnetic vertex that they 
acquire in the presence of matter. We calculate the contribution to the photon self-energy in a 
, dense plasma, due to the presence of a gas of charged particles, or neutrinos, that moves as a whole 

• relative to the plasma. General formulas for the transverse and longitudinal components of the 

' photon polarization tensor are obtained in terms of the momentum distribution functions of the 

04 ■ particles in the medium, and explicit results are given for various limiting cases of practical interest. 

f—j ' The formulas are used to study the electromagnetic properties of a plasma that contains a beam of 

I neutrinos. The transverse and longitudinal photon dispersion relations are studied in some detail. 

Our results do not support the idea that neutrino streaming instabilities can develop in such a 
system. We also indicate how the phenomenon of optical activity of the neutrino gas is modified 
due to the velocity of the neutrino background relative to the plasma. The general approach and 
results can be adapted to similar problems involving relativistic plasmas and high-temperature gauge 
theories in other environments. 



I. INTRODUCTION AND CONCLUSIONS 



> 

o 

, From a modern point of view, the methods of finite temperature field theory (FTFT) M provide a natural setting to 

■ treat the problems related to the propagation of photons and neutrinos through a dense medium. This view has been 
[ largely stimulated by the work of Weldon (^-Q], who emphasized the convenience of carrying out covariant, real-time 

(~| ■ calculations in this kind of problem. The work of Weldon demonstrated that the real-time formulation of FTFT is 
rS I well suited to the study of systems involving gauge fields and/or chiral fermions at finite temperature, which can be 
JL , extended in an efficient and transparent way to realistic situations involving, for example, photons and/or neutrinos 
1^ in a matter background. 

, The electromametic properties of neutrinos in a medium, besides their intrinsic interest, are relevant in many phys- 
l!" ■ ical applications [|6|. For example, the induced electromagnetic couplings of a neutrino propagating in a background 
. !^ I of electrons and nucleons is responsible for the plasmon decay process 7 — > i/P in stars, and modify the MSW resonant 

■ condition in the presence of an external magnetic field @-|l^] . A neutrino gas also exhibits the phenomemon of optical 
5-H I activity as a result of the chiral nature of the neutrino interactions 

^ The covariance in this type of calculation is implemented by introducing the velocity four-vector of the medium, 

in terms of which the thermal propagators are written in a manifestly covariant form. Therefore, covariance is 
maintained, but quantities such as the photon self-energy or the neutrino electromagnetic form factors depend not 
just on the kinematical momentum variables of the problem, but also on u'^. In practice the vector is in the end set 
to (1,0), which is equivalent to having carried out the calculation from the start with respect to a frame of reference 
in which the medium is at rest. This is usually the relevant situation. Therefore, while generally useful, the covariant 
nature of these methods has not been of particular importance in the applications mentioned. 

However, there is a class of problem in which setting m'' — (1,0) is not possible. These are problems that involve 
one stream of particles (which we can think of as a moving medium) fiowing through a background medium, which 
we can take to be stationary. If we denote by it^ the velocity four- vector of the stationary medium, and by u' the 
corresponding one for the moving background, then we can set u'^ = (IjO), but we cannot take both to be (1,0) 
simultaneously. Thus, for example, if we were to calculate the self-energy of the photon propagating through such 
media, it will depend on the momentum and velocity four- vector u'^ as usual, and in addition on u' . This additional 
dependence can have consequences that are as important as the effects of the stationary background itself. 

For example, it is well known that in a plasma in which a bunch of electrons move, as a whole, relative to a plasma at 
rest, besides the usual dispersion relation of the longitudinal photon mode, another branch appears whose dispersion 
relation depends on the velocity of the beam. Under some conditions, the sign of the imaginary part of this dispersion 
relation is such that the corresponding amplitude grows exponentially, which signifies an instability of the system 
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against the excitation of those modes. This kind of system is famihar in plasma physics research, and examples of 
them are discussed in textbooks on the subject 13 1^. 



Recently ]15( | , it has been sug gest ed that a similar kind of streaming stability might be driven by a flow of neutrinos 
through a matter background ||l^ , |l7|] . Because the neutrino acquires an effective electromagnetic coupling as it 
traverses a medium p^ , |l^ ,[7|,p^ , the propagation of a photon in a medium that contains a drifting neutrino background 
may be affected in a way similar to the case mentioned above. As argued in Ref. |p^ , such effects can appear under 
the conditions of realistic situations such as those in a supernova explosion, gamma ray bursts, or the Early Universe. 

Similarly, other neutrino processes that have been studied previously, such as those mentioned above, may be 
modified in important ways if the neutrino gas is moving as a whole relative to the matter background. 

Motivated by all these considerations, in this work we calculate the neutrino contribution to the photon self-energy 
in a medium in which the neutrino gas moves as a whole relative to a matter background which we take to be at 
rest. The calculation is based on the application of FTFT to this problem in the manner that has been suggested 
above. The implicit assumption is that, in the rest frame of the stream, the neutrino background is characterized 
by a momentum distribution function that is parametrized in the usual way. Although our focus is the case in 
which the neutrino background constitutes the stream, largely motivated by the potential applications that have been 
mentioned, the calculation and the formulas for the photon self-energy are presented in such a way that they can 
be adapted to other cases as well. Therefore, they complement the existing calculations of the photon self-energy in 
which all the particles form a common background with a unique velocity four-vector. The results for the photon 
self-energy can be equivalently interpreted in terms of the dielectric constant of the system, and in that way we 
show that the well known textbooks results for the stream stabilities are reproduced when the appropriate limits are 
taken. On the other hand, the results we obtain are valid for general conditions (whether they are relativistic and/or 
degenerate or the converse) of the gases that form the plasma at rest as well as the stream, hold for general values of 
the velocity of the stream, and are valid also for general values of the photon momentum and not necessarily for some 
particular limit. Therefore, they are useful also in the study of similar processes that may occur in other contexts, 
such as high-temperature QCD, heavy ion collisions or other similar environments in which the methods of FTFT 
are applicable. 

In Section 0, we give the general one-loop formulas for the generic contribution of a moving fcrmion background 
to the photon self-energy. The contribution from any given fcrmion can be written in terms of a few independent 
functions, which are expressed as integrals over the momentum distribution functions of the fermion. Explicit formulas 
are given for various limiting cases of physical interest, which also serve to show how some of the results derived in 
textbooks for simple cases are reproduced in the appropriate limit. 

The case of the system that is composed of a matter background made of electrons and nucleons (and possibly 
their antiparticles), and a neutrino gas that propagates as a whole relative to the matter background, is considered in 



Section [II . We begin by reviewing the one- loop formulas for the effective electromagnetic vertex of the neutrino in a 
matter background, in the way that will be used in the calculation of the photon self-energy. The neutrino background 
contribution to the photon self-energy is then determined. It depends on the momentum distribution functions of 
the matter particles and well as those for the neutrinos. As an application of the formulas obtained, the dispersion 
relation for the longitudinal modes is considered, with attention to the possible effect of the neutrino contribution 
to the instability of the system. In that context, our results do not indicate the existence of unstable modes, and 
therefore we do not find support for the idea that stream instabilities due to the presence of the neutrino background 
can develop in such systems. 



In Section III we also consider the dispersion relations for the transverse modes. The chiral nature of the neutrino 
interactions gives rise to the phenomenon of optical activity, which had been studied earlier . Here we show how 

the results of Refs. pT| , p^ are modified when the neutrino gas is moving relative to the matter background. The main 
effect is that the dispersion relations for the two circularly polarized modes are not isotropic. As a consequence, the 
frequency splitting between them, which is the measure of the rotation of the plane of polarization, depends on the 
direction of propagation of the mode relative to the velocity of the neutrino gas. In particular, under the appropriate 
conditions, the frequency difference is the opposite to what is found if the neutrino gas is not moving relative to the 
matter background. 



Section IV contains our outlook, where other possible effects and applications arc also mentioned. 



II. PHOTON SELF-ENERGY IN A FERMION BACKGROUND 

We will consider a medium that consists of a gas of nucleons, electrons, neutrinos and their antiparticles. Each 
fermion gas gives a contribution to the elements of the 2x2 photon self-energy matrix, that are determined by 
calculating the diagram shown in Fig. [I]. 
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p + q 

inll^^ (q) = q \rJv!A(a^^^ q 
P 

FIG. 1. Diagram for the contribution to the photon self-energy matrix from a generic fermion. For a charged fermion, the 
electromagnetic coupling is given by the tree-level terms in the Lagrangian while for the neutrino it is the one-loop vertex 
function induced by the matter background. 

In particular, the contribution to the tth^j/ element from each fermion / in the loop is 

^-Ki. = (-l)(-O^Tr J ^4;")(,)^4/)^(p-,,)47)(-,)^4/\(p), (2.1) 
where j'f^^q) is the electromagnetic current of the fermion. It is defined in such a way that the on-shcU matrix 



element of the electromagnetic current operator jjf"^"^ is given by 

(/(p')I^('^"'(o)mI/(p)> = u{p')jf;;\q)u{p) , 



(2.2) 



where q = p — p' and u{p) is a Dirac spinor. For the electron it is simply 67^, for the nuclcons it must in principle 
include the magnetic moment term, and for the neutrino w e mu st use the effective electromagnetic neutrino vertex in 
the medium. The fermion propagator that appears in Eq. (2.1) is given by 



4n(p) 



where 



{i> + mf) 



2m5{p^ 



'f)Vf{P' 
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(/)- 



+ 1 



+ 



0{-p- 
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+"/ + 1 



(2.3) 



(2.4) 



with being the inverse temperature and a/ the fermion chemical potential. The vector is the velocity four- 

vector of the fermion gas as a whole, so that u^^^'^ = (1, 0) if the fermion background is at rest. In Eq. ( ^.4|) we are 
allowing for the possibility that the different fermion gases of the background may be at different temperatures and, 
most importantly for our purposes later, that each gas has a velocity four-vector that is not necessarily the same for 
all of them. The implicit assumption here is that, in the rest frame of each fermion background, the corresponding 
particles have an isotropic thermal distribution characterized by a temperature and chemical potential 1//3/ and a / . 
The dispersion relations of the propagating photon modes are obtained by solving the equation 



where 



/ 



(2.5) 



(2.6) 



and we have denoted by ir^fj the background-dependent term of Eq. (^Tj). In the rest of this paper we will focus only 
on the real part of the dispersion relations, but similar considerations could be used to calculate the imaginary part 
as well. 

In order to calculate Re7r^1?^\ and thus determine the dispersion relations, we must know the composition of the 
medium and the formulas for the electromagnetic current that must be substituted in Eq. (2.1). To proceed, we 
consider the various cases separately. 
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A. Matter background 



We consider first an isotropic medium composed of nucleon and electron gases, with a common velocity four-vector 



u^. The most general form of the physical self-energy function in this case, which we denote by tt^™"* is |2|,|1 1| 



(2.7) 



where 



R^uiq, u) = -2 Qp.u{q, u) 



P^^u{q■, u) = —e^,l,a|3q'^u'^ , 



with 



Although we have not indicated it explicitly, in general, tt,^^ p are functions of the scalar variables 

~ q ■ u 



{m 



(2.8) 
(2.9) 

(2.10) 



which have the interpretation of being the photon energy and momentum in the rest frame of the medium. 

A detailed calculation of the photon self-energy in such a medium was carried out in Ref. [ pO| . For our present 
purposes, it is useful to summarize those results as follows. The nucleon magnetic moment term contribution is not 
important for practical purposes. Therefore, we use here jj^™'' = 6/7/^; so that the neutron contribution is being 

neglected. Substituting in Eq. (2T) the formula for the contribution from each fermion in the loop can be 

expressed in the form 



Re 7r(^{,) = -46^ 



Rf,i,{q,u) + ^^^^^Qp.i'iq, u) 



(2.11) 



with the coefficients Af and Bf defined as 



Bf{n,Q)^ 

In these formulas. 



[ff{p-u) + fjip-u) 



2toj — 2p- q 



+ iq^ -q) 



q^ + 2p ■ q 

2{p ■ uf + 2{p ■ u){q ■ u) - p ■ q 
q'^ + 2p ■ q 



+ iq^ -q) 



p 2 + mi , 



and /j J denote the particle and antiparticle number density distributions 



(2.12) 



(2.13) 



(2.14) 



with the minus(plus) sign holding for the particles (antiparticles), respectively. Comparing Eqs. (2.7) and (2.11) we 
can identify the contribution of any fermion to the real part of the transverse and longitudinal components of the 
self-energy, and therefore obtain 

Re4'")=.-2e^^fA;(f^,Q) + ^i?/(f7,Q)) , 



Re4"'=4e2 5:|,%(f^,Q), 



(2.15) 



where the relation = —Q^/q^ has been used. 
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B. Matter background and a stream of charged particles 



We now consider a medium that contains, in addition to the background as has been considered above, another 
gas of particles with a velocity four- vector u'^. We will refer to them as the matter background and the stream, 
respectively, and we assume that the latter consists of only one specie of fermions /' with an electromagnetic coupling 
jj.?™-* = e/'7^. The fermion /' could be, for example, the electron or any other charged particle. We will denote by 
C/'° and U' the components of u'^ in the rest frame of the matter background so that, in that frame. 



u^ = (l,0), m'^ = (C/'°,C/') 



The contribution from /' to the photon-self-energy is given by a formula analogous to Eq. (2.11), 
where 



(2.16) 



(2.17) 



Re TT, 



(/') 



(2.18) 



In Eq. (2.18) the functions Afi and Bfi are given by formulas ana logo us to Eq. (2.12), but with the replacement 



u'^, and we have used u = —Q /q where, similarly to Eq. (|2.9|) , 



(2.19) 



In analogy with Eq. (2.10), the variables O', Q' are defined by 

Q'^^n'^-q\ 
and they are expressed in terms of and Q by the relations 

= u'°n-u' ■ Q 



(2.20) 



Q' = \ {u'"n-u' ■ Q] -n'^ + Q 



The total photon self-energy is given by 

= 4'"^i?^.(g,7/) -f 4"^Q^.(q,u) +4'^i?^.(g,7/') +4'^Qm-('Z,"')- 



(2.21) 



(2.22) 



Eq . (2.22) can be written in a convenient form by the following procedure. From the definition of R^i, given in Eq. 
(2.8) we have 



Rfj.viq, u') = Rtj.,yiq, u) + Qfj,u{q, u) - Qfj.u{q, u') . 



We now define the vectors 



with 



^_ R^^''{q,u)u', 



= -^P^"(g,w)el,. 



(2.23) 



(2.24) 



{u-u'f ^.^ 



(2.25) 
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which can be expressed in the form 



It is easy to verify that e'^ ^ are mutually orthogonal and satisfy 



ei,2 • g = ei,2 • M = , 



p2 



-1. 



(2.26) 



(2.27) 



Thus, together with u^, they form a complete set transverse to g^, and therefore i t is p ossible to express u'^ in terms 
of them. The desired relation, which follows from substituting Eq. ( ^.^ ) into Eq. ( 2.24 ), is 



which substituting in the definitions given in Eq. (^^) yields the convenient formulas 



(2.28) 



g/^2 ^A^-v^'-/ ■ g,~2 

On the other hand, as shown in Rcf. R^^ can be decomposed in the form 

Rtj,u{q,u) = -(e^ei 



=2^2^ 



(2.29) 



(2.30) 



In this way, using Eqs. ( 2.23| ) and ( 2.29|) in Eq. ( 2.22 ) together with the decomposition given in Eq. ( 2.30| ), we finally 
arrive at 



(m) 



m (if 



£2^62 



1/ (tt, 



(m) 
T 



(mj_mO / (/') 



(2.31) 



Eq. (2.31), besides unfolding the main structure of the modes in a particularly simple way, is a useful formula that 
allows us to obtain the dispersion relation of the modes under a variety of conditions. In the absence of the stream, 
the solutions consist of one longitudinal mode with polarization vector 63 oc m'^, and two degenerate transverse modes 
with polarization vectors e^2 that satisfy Q^vS.i2 = 0. Their dispersion relations are determined by solving the 

equations (f' — Re 7rj^"^ for the longitudinal and transverse modes, respectively. The presence of the stream breaks 
the degeneracy of the transverse modes, and in general causes a mixing between the m wit h the longitudinal mode. 
In those cases in which it is permissible to treat the mixing term (the last term in Eq. (2.31)) as a perturbation (e.g., 
the number density in the stream is sufficiently smaller than those in the matter background), the dispersion relations 
are obtained approximately by solving 



q' 


= 4")- 


h4') + 


q' 






q' 


- 4"^ - 


,4')- 



g/2 

Niq" 
Q 



M') 



.5.") 



7^(4'"' 



J/') 



(2.32) 



with corresponding polarization vectors ei.2 and 63 oc m, respectively. In Eq. ( 2.32| ) we have used the relation 
= —Q^/q^ and the corresponding one for u'"^. If the mixing term is not sufficiently small so that the higher order 
terms are important, then the full 2x2 problem in the ei — u plane must be considered which, although tedious, is 
straightforward. 

In the equations Eq. ( ^.32| ), it is understood that the variables il', Q' are to be expressed in terms of f2, Q by means 
of the relations given in Eq. (2.21). They thus become implicit equations for f2, Q, whose solutions determine the 
dispersion relations il(Q) of the various modes. 
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C. Discussion 



For illustrative purposes, we consider the specific case of a stream of electrons and a matter background that consists 
of an electron gas and a non-relativistic proton gas. We borrow from Ref. jl^ [see Eqs. (A5) and (A9)] the following 
results 



1 f d'V Q-\7v{ff{S)+fj{S)) 



(27r)3 



S-p vvVv{ff{£) + fj{£)) 



{2n) 



(2.33) 



V-p is the gradient operator with respect to the momentum V and v-p — V /£. As shown 



where £ = 

there, they are valid for values of q such that 



ql{£) « 1 : 



(2.34) 



where {£) stands for a typical average value of the energy of the fermions in the gas. For distribution functions that 
depend on V only through £, we can replace V-p V'p-§g in Eq. (2.33). Several useful formulas follow from Eq. (2.35) 
in particular cases. For example, if the fermions are relativistic. 



Of 



Of 



1 n 

2Q 



n-Q 



(2.35) 



Eq. (2.35) holds for a degenerate or non-degenerate gas. For the non-relativistic and non-degenerate case we use 



Bf{n, Q) 



-3woy - 

f]2 



Of 



which are valid if, in addition to Eq. ( 2.34 ), 



n > vfQ, 



(2.36) 



(2.37) 



where Vf = 1/ yj iijvrij is a typical value of the velocity of the fermions in the gas. The quantity ojqj, which is related 
to the plasma frequency in the gas, is given by 



^0/ 



:(//(£) + //-(£)) 



3f2 



(27r)32f 

In Eqs. (IHI) and (|3^) we have used its form in the relativistic (ER) limit and non-relativistic (NR) limits, 

^ Jo°° dVViffiV) + fj{V)) (ER) 



(NR), 



(2.38) 



(2.39) 



where Uf is the fermion number density in the frame in which the background is at rest. The corresponding formulas 
for the non-relativistic and degenerate case are given in Ref. |po| . 

Under most circumstances, the protons make a negligible contribution to the dispersion relations. The conditions 
under which those terms can be important are given in Ref. po|| . Here we do not include those special cases and 
therefore we will not consider the protons further. Lastly , we assume that the stream is not moving too fast as a 
whole, so that the term that is proportional to A^i in Eq. ( 2.32| ) can be neglected, since according to Eq. ( 2.26 ) it is 
or the order of the velocity squared of the stream. 

With these assumptions, the dispersion relations become 



q^ = 4e2 



Ae'in\Q') + -^Be>in',Q') 



^^B,{n, Q) + -^BAn',Q') 



(2.40) 



for the transverse and longitudinal modes, respectively. We now consider several cases separately. 
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1. Non-relativistic matter and stream electrons 



For the electrons in the matter background we use Eq. ( 2.3(^ ). Similarly, for the stream 

n'2, ,2 



n'2, ,2 



f^'2 



(2.41) 



which, as we will see, are suitable for finding the long wavelength limit of the dispersion relations. Eq. ( ^.41 ) can be 
used if the solution is such that 



n' > Ve' Q! . 



(2.42) 



The conditions under which the solution thus found is valid can be ascertained afterwards. The formula for uj^^, is 
the same expression given in Eq. ( ^.38 ), but with replacement ff j{£) fe\e'{£)^ where 

1 



(2.43) 



As we have indicated previously, the implicit assumption that is being made here is that the electrons that compose 
the stream have, in the rest frame of the stream, an isotropic thermal distribution characterized by a temperature 
and chemical l//?e' and Ue' ■ 

Let us consider the dispersion relation for the longitudinal mode. From Eqs. ( ^.15 ) and ( ^.18 ) this yields 



1 = 4e^ 



r!2 



From Eq. (|2^ 

using U'^ ~ 1, and therefore the dispersion relation is 



0. 



(2.44) 



(2.45) 



(2.46) 



The salient feature of Eq. ( 2.46| ) is that, besides the usual solution 0|(Q ^ 0) ~ ie'^ujQ^, there is another one with 
~ U' ■ Q. The standard way to find this second solution is to substitute 



n:^U' .Q + 6l 

in Eq. ( ^.46| ) and determine approximately by taking it to be a small quantity. In this fashion, we find 

\2eu;oe'U' ■ Q\ 



± 



'((7'- Q)2-4e2u;2^ 

which shows the well-known instability of this system. For values oi U' ■ Q such that 

0<|C7'-Q|<2|e|cj0e, 



(2.47) 



(2.48) 



(2.49) 



the dispersion relation has a solution with a positive imaginary part, which signals that the system is unstable against 
oscillations with those values of U' ■ Q. The condition that 6l be small relative toU' -Q is satisfied for sufficiently small 
values of woe'/^Oe- On the other hand notice that, for this solution, il' = 5l, and Q' = y/ Q' — fl"^ + Q'^ ~ Q. The 
conditions given in Eqs. (2.34) and (2.42) are then equivalent to \U' ■ Q\ ^ VeQ and S ^ Ve'Q which, for sufficiently 
small values of the thermal velocities, are satisfied as well. 

Turni ng no w the atten tion to the tran sverse dispersion relation, we substitute the formulas for Ae^gi and i3e,e' given 
in Eqs. ( p^ ) and ( ^.4lD into Eq. ( ^.40| ). This yields simply 



9' = ^eW^^ + 4eWo^, 



(2.50) 
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which shows that in this case the presence of the stream perturbs somewhat the transverse dispersion relation by 
shifting the v alue of the plasma frequency, but it does not produce a significant effect otherwise. 

Eqs. (2.44) and (2.5C) reproduce the well-known results found in textbooks ||2^, which are derived by kinetic theory 
or sim ilar s emi-c lassical methods. However, the results that we have obtained, and which are summarized in Eqs. 
(2.31) and (2.32), go fa rther. Toge ther with the expressions for the self-energy functions in terms of the coefficients 
Af and Bf [Eqs. ( 2.15 ) and ( 2.18| )] they allow us to study systems under a wider variet y of conditions, including 
those for which the semi-classical approaches, and the simple formula given in Eq. ( ^.44 ) in particular, are not be 
applicable. 



2. Non-relativistic matter electrons and relativistic stream electrons 



For the electrons in the stream we must use in this case 



Oe' 



(2.51) 



while the matter electron formulas are the same as the previous ones. The dispersion relations are then determined 
by 



C!^-4e^<=4e^c.o^,,/L 
for the longitudinal and transverse modes, respectively, where we have defined 



(2.52) 
(2.53) 



2Q' 



■In 



2 r ■ 

Let us consider the longitudinal dispersion relation. 



1 - ^^In 
2Q' 



1 



n' + Q' 



n'-Q' 



(2.54) 
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U' = 0.2,61 = 45' 
U' = 0.5,6 = 



2 4 6 8 10 




n/Q 

FIG. 2. The functions /l,t defined in the text, are plotted as functions o{ Q/Q, for some representative values of the velocity 
of the stream U' and the angle between Q and U' . For Q/Q — 1 the function /_l becomes infinite, while fx ~ 3/2. 

Wi th the help of Fig. || it is easy to see that, besides the usual solution il^ — 4e^cjoe, any other solution to Eq. 
(2.52) must have SI « Q (which implies that fi' « Q' and the function can be large) so that the stream term can 
compete with the matter term in that equation. Substituting 



in Eq. (|^) we find 



n = Q + 5L 



2 2 



(2.55) 



(2.56) 



for > 4e^(jjQg. In contrast to the case considered in Section II C 1, there is no sign that a stream instability may 
develop in the present one. 

For the transverse dispersion relations the situation is different. The function Jt is not larger than a number of 
order 1, as shown in Fig. ^, so that the stream contribution in Eq. (2.52) only introduces a perturbation in the usual 
solution. 

In summary, when the stream consists of a relativistic electron gas, there is no sign that any stream instabilities 
may develop. This result will be a useful reference point when we consider in Section III the case in which the stream 
consists of neutrinos. 



III. PHOTON SELF-ENERGY IN A NEUTRINO STREAM 



In this section we consider a system composed of a matter background as in Section [I B , and a neutrino stream 
with a velocity four- vector w'**. Our immediate task is to determine the neutrino stream contribution to the photon 
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self-energy, for which we must calculate a diagram similar to the one in Fig. but with a neutrino as the fermion in 
the loop. Denoting the effective electromagnetic vertex of the neutrino in the matter background by T^ii^\q), then 



(3.1) 



where the neutrino propagator S^pl^ is is given by Eq. (|]^) (with rriu = 0). The neutrino effective vertex can be 
expressed in the form 



where L = i(l — 75) as usual, and T^^ can be decomposed as 



(3.2) 



(3.3) 

A detailed calculation of the various terms in Eq. ( |3.3| ) was carried out in Rcf. ||10|]. For our purposes here, we can 
summarize the main results obtained there as follows. 

For practical purposes, the contribution to Tt,l due to the anomalous magnetic moment couplings of the nucleons 
in the background is negligible, so that 



while 



Tt = 2V2\e\GFa, (^A,{n, Q) - ^^^^^^ + , 
Tp = - 4V2GFfopQ(|e| + 2mpKp)Cp{n, Q) - 8m„K„\/2Gi.fo„QC„(17, Q) . 



(3.4) 



(3.5) 



In these formulas 



Kp = 1.79 



2mr 



(3.6) 



are the anomalous magnetic moment of the nucleo ns, th e coefficients ay and bf are the neutral current couplings of 
the fermion /, while Ap and Bp are defined in Eq. ( 2.12| ) and 



u ■ p 



iff ~ ff) 



1 



q^ + 2p ■ q 



+ {q^ -q) 



(3.7) 



The electron terms Tj^' were calculated in Ref . 0| , and are given by 



4^^^2V2eGp (^A.(f^,Q)-M^^| 
^.(^)=4^/2eG.M^ 



fle + 1 (lye) 
a-e (t'M.r 



T^"^ = -4V2eGFQCe{n,Q) 



fle + 1 (j/g) 
be - I {Ve) 



(3.8) 



Substituting Eq. (3.3) in Eq. ( p.l| ), wc then obtain for the neutrino stream contribution to the photon self-energy 

Rc^^-^,) = -2T^„(g)T,^(-g)J"/' (3.9) 
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where 



TQ/3 ^ 



1 



~ 2p ■ q 
1 



+ {q^ -q) 



q^ — 2p ■ q q^ — 2p ■ q 



(3.10) 



The transversality and symmetry properties of J"'^ imply that it is expressible in t erms of the tensors R"^ {q,u'), 
Q°'^{q, u') and P°'^{q, u'), with coefficients that can be determined by projecting Eq. ( 3.10 ) along these tensors. Thus 
we find 



+ ^^4^0"'^(g, u') + Q'CA^', Q')P"^iq, u') , 



(3.11) 



with the coefficients A^i^l' , Q!), B^{n', Q') and C^(fi', Q') defined in Eqs. (^.12[ ) and (U). In Eq. (|1|) the tensors 
R^p(qyu!), Q^i,{q,u') and Pfj_^{q,u') are eliminated using Eqs. (2.23) and (2.29), and R^^(q,u) is decomposed as in 
Eq. (|2.30[). In this way, 



<^ a/: 



2 I u'2 j 2m'2 



2w'2 



-A, 



Qaf3{q,u) + CyQ 



Paf}{q,u) 



+ 



2w'2t22 



3B, 



(eiQU/3 



+ Mae 



1/3) + 



(£(0.62/3 - e2aU/3) . 



(3.12) 



By substituting Eq. (3.12) in Eq. ( |3.9[ ) we finally obtain the formula for the neutrino contribution which, when it is 

added to tt^™'' to obtain the total photon self-energy, is the starting point to determine the photon dispersion relations. 
However, with all its generality, the formula is not particularly useful and therefore we consider below some specific 
situations of potential interest. 



A. Longitudinal dispersion relation 



As already seen in Section |1I B| , in general the effects of the stream break the degeneracy of the transverse modes 
and also mixes them with the longitudinal one. When the latter effect is not too large, the longitudinal dispersion 
relation is obtained approximately by solving the equation 



(3.13) 



where 



(3.14) 



Using the relation u^Tf^a ~ T^Ua, together with 



-Rf,i,{q,u') = --^72 



}tiu{q,u') 



(3.15) 



we obtain from Eq. (3.9) 



12 



ReTT, 



Q 



12 



-A, 



(3.16) 



where, to simplify the notation, we have omitted the arguments 57' and Q' in the coefficients Ai, and B^. The 
dispersion relation obtained by substituting Eqs. ( 2.15| ) and (3.16) in Eq. (3.13) is the same as the corresponding one 



for a stream of charged particles, with an effective charge 



V2 



1. Long wavelength limit 



We consider the case analogous to the one discussed in Section II C , namely, a matter background made of a non- 
relativistic electron gas and a non-relativistic nucle on ga s, and the neutrino stream. As in the case mentioned, the 
photon momentum is assumed to be such that Eq. (2.34) h olds. For the electrons in the matter background we then 
use th e formulas for Ae{^^ Q) and Be{^, Q) given in Eq. ( [2.36| ), which are valid for fl ^ v^Q a.s indicated in Eq. 
(2.37), while the analogous proton terms are negligible. From Eqs. (3.4) and (3.8), this yields 



(3.17) 



where 



To = 4V2G 




(3.18) 



The neutrinos are, for all practical purposes, ultrarelativistic. For them, we use the formulas analogous to those in 
Eq. (2.35), namely 



B,in',Q' 



-3 a;, 



2Q' 



In 



n' + Q' 



(3.19) 



with 



1 

6^ 



dVViUiV) + MV)) ., 



(3.20) 



where 



(3.21) 



are the neutrino and antineutrino momentum distributions, in the rest frame of the stream. The longitudinal dispersion 
relation obtained from Eq. (3.13) is 



2T^ 



(3.22) 



In this equation, Q' and Q' are expressed in ter ms of Q and Q by means of Q' = ^fi' ^ — with ft' = U'^^l — U'- Q, 
as indicated by Eq. ( ^.20 ). The solutions of Eq. ( 3.22| ) determine the dispersion relation VIl{Q,) in the long wavelength 
limit and are valid for f2 3> VeQ- 



B. Neutrino driven stream instabilities 



Besides the usual solution Vl\ 



4e Woe: Eq. (3.22) can have an additional solution if the neutrino term is sufficiently 
large that it can compete with the electron term. To determine whether this can happen, cons ider the specific situation 
in which the velocity of the neutrino stream is not too large, so that the term in Eq. ( ^.22| ) proportional to Ni can 
be neglected [see Eq. ( ^.26 )]. Using the formula for By given in Eq. ( 3.19| ), the longitudinal dispersion relation then 
becomes 
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(3.23) 



where is the same function defined in Eg. (2.54). For values of w Q the function Jl becomes large, but with the 
factor (1 — Q^/ri^) its contribution in Eq. (3.23) is negligible in that region. On the other hand, for ^ 0, 



so that 



(3.24) 



(3.25) 



in this limit, instead of Eq. (3.17). Therefore, the neutrino contribution is not large in the limit il either. 

We therefore conclude that the neutrino contribution produces a small correction to the usual dispersion relation 
but does not introduce any additional branch. In particular, there are no stream-induced instabilities in this system. 
This conclusion is in sharp contradiction with the result obtained in Ref. where it was found that, in the same 
system, the dispersion relation indicates the appearance of neutrino-driven stream instabilities. 

To understand the origin of this discrepancy it is useful to consider 



(3.26) 



for the momentum distribution function of the neutrinos, which is of the form used in Ref. |15|. Using it in Eq. ( 2.33| ) 
to calculate B,, results in 



2£ 



Q^-{ Q- U' f 



which, when substituted in Eq. (3.22), yields the longitudinal dispersion relation 

>2 o2\2 g2gij^2, 



f7"(f]' -4eXe) 



To{nu + ftp) 



£ 



(17- Q cos 61)2 ' 



(3.27) 



(3.28) 



where cos Q 
form 



with 



Q • V' and we have neglected the term proportional to A^i, as before. Eq. ( 3.281 ) has a solution of the 



= Qcos6' + 5 



Av-)i _ Tpin^ + np) 



Q2 sin^ 6 



£ 



C0s2 0(Q2 cos^ 



0-4e2a;2j' 



(3.29) 



(3.30) 



which exhibits an instability for cos^ ( 

Thus, while we are able to reproduce qualitatively the result of Ref. [|15[ in this way, we must note that it is based 
on an inconsistent application of the long wavelength formulas given in Eq. (2.33). As explained in detail in Re f. P0[ , 
those formulas are obtained by expanding the integrands in powers of q/£ in the one-loop formulas given in Eq. ( p. 121 ), 
and retaining only those terms that are dominant in the limit q/ £ 0. This requires, among other conditions, that 
the momentum distribution functions be such that its derivatives do not introduce any singularities in the integrands. 
The results derived in this way are equivalent to those obta ined by semiclasical methods based on kinetic theory 
or similar approaches. However, the form given in Eq. (3.26) does not satisfy the required conditions and therefore 
neither the long wavelength approximation of the one-loop formulas, nor the semiclassical formulas, are applicable in 
that case. 



Leaving aside the question of whether or not a distribution function such as that given in Eq. (3.26) is realistic 
in any particular physical cont ext, i n order to use it the coefficients A^^B^ must be calculated with the complete 
one- loop formulas given in Eq. (2.12). Following this procedure we obtain 



B„ 



2£{n^ 



l\2 



4£2(J7 _ Q . U'Y - (fi2 - Q2 



(3.31) 
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instead of Eq. ( |3.27 ), so that the longitudinal dispersion relation becomes 



Q2 sin^ 



(3.32) 



If we neglect here the term in the denominator, we recover Eq. (3.28). However, when that term is taken into 
account , the neutrino contribution does not become large for « Q ■ U' , and therefore a solution of the form given 
in Eq. ( ^.2S| ) does not exist. 



C. Transverse dispersion relation 

The dispersion relations for the transverse modes are given approximately by solving the equation 



/ 2 (™) \ r> t \ 'y'^j 



where 



(3.33) 



= -2 [TTiq)R^x{q, u) + Tp{q)P^x{q, u)] [TT{q)R,p{q, u) - Tp{q)P,p{q, u)] J^P 



(3.34) 



i s th e transverse projection of the neutrino contribut ion t o the self-energy, and in the second line we have used Eq. 
(3.9). For J"^ we will use the formula given in Eq. (3.11) and consider the case in which the terms with the factor 
A^i can be neglected, as we did in Section [II B. Therefore, remembering Eq. (2.30), we will substitute in Eq. (3.34) 



\ M--|^)^a/3(9:") + C^-2(^)^"/3(g,u) 



It is now useful to introduce the linear combinations 



which satisfy 



and have the representation 



where 



= \ {Rali{,q-,U) ± Pai3{q,u)) 



D(s)a/3D(s') _ c ca 



R 



(±) 

'al3 



(3.35) 

(3.36) 

(3.37) 
(3.38) 
(3.39) 



Writing Raf3 and Pap in terms of R^^'' and substituting the resulting formulas in Eq. (3.34), with the help of Eq. 
( 3.37 ) we obtain 



where 



7r(±) = -2 (Tt ± Tp) 



From Eq. (3.33), the dispersion relations for the transverse modes are then 



(3.40) 
(3.41) 

(3.42) 



with the corresponding polarization vectors Cq^'. 
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D. Optical Activity of the neutrino gas 



Eq. ( ^.41 ) exhibits the phenomenon of optical activity of the neutrino gas, in which the two circularly polarized 
photon modes travel with different dispersion relations as a result of the chiral interactions of the neutrino 
Notice that for this occur, Tp and/or must be non-zero. This requires that the chemical potentials in the matter 
background be such that, for some particle species, a/ 7^ 0, or that ^ 0. In the latter case however, there is an 
additional contribution to the photon self-energy that arises from the set of diagrams that were calculated in detail 
in Ref. ]T^ . Those diagrams are not included in Fig. |l| and their result is an additional contribution to the photon 
self-energy that must be taken into account in Eq . (3.42 ). The result of the calculation of Ref. jl^ is taken into 
account by including in the right-hand side of Eq. ( |3.9D the term 



(3.43) 



where, in the notation of the present paper, 



(3.44) 



with 



(2^ 



(3.45) 



The result quoted in Eq. (3.44) is valid for values oi q < me- When this term is included in Eq. (3.34), the net effect 
is that the right-hand side of Eq. ( 3.42 ) has the additional the term 



Q \ u-u' 



-1 



(3.46) 



If the background contains an equal number of neutrinos and antineutrinos, then Hp'' a s wel l as are zero. In such 
a case, the optical activity of the neutrino gas is due to a non-zero value of Tp in Eq. (3.41), which in turn depends 
on the difference between the particle and antiparticle number densities in the matter background. 



1. Long wavelength limit 



TT^^-* can be evaluated explicitly by considering specific situations. As an example we consider once more a matter 
background that consists of non-rclativistic electron proton gases, with the photon momentum satisfying Eq. (2.34) 
and Q 3> VfQ. The proton contribution to Tt is negligible, and using Eq. ( p. 33 ) for Ae and Be, 



Tt — —To , 



(3.47) 



with To given in Eq. (|t|). On the other hand, Tp is given by Eq. ( |3.5[ ) where, in the momentum regime that we are 
considering, 



1 f dPV fff-ff 



2 J (27r)32f V £ 



2V^ 
3f2" 



(3.48) 



as shown in Ref. For the electron and proton gases we use the non-relativistic limit of this, C/ = 

which implies that the proton term is negligible and therefore 

Q 



''Of 



/2m/, 



where 



Tf 



T' = 4V2G 



2me 



T' 




(3.49) 



(3.50) 
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For the neutrino gas, the relativistic hmit of Eq. (3.48) yields 

1 



247r2 J, 



dV{MV)-MV))., 



(3.51) 



while Ai, and are given in Eq. ( 3.19| ) 

With the help of these f 
dispersion relation becomes 



With the help of these formulas and remembering that ReTr^"-* — Ae^u^^ in the case we are considering, the 



2 A 2 2 



It 



n 2 r 



± 



fu-u' 

Q! 



(3.52) 



where we have included the Hp term as indicated in Eq. ( 3.46| ). 

Let us consider first the situation in which k. /p, so that Hp'' and Ci, can be neglected in Eq. ( 3.52| ). The 
solutions for the two modes are then given by 



17± 



Q2 + 4e2w^g =F 



which is of the same form as that given in Eq. (4.22) of Ref. |ia|, if we make the correspondence 



(3.53) 



(3.54) 



there. For the situations of potential interest analyzed in Ref. ||T^, the effects of the dispersion relations given in 
Eq. ( 3.53| ) are smaller than those found in that reference by a factor of about Gpuj^^ ~ Gpne/me^ and whence are 
unimportant for all practical purposes. 

Therefore, retaining only the term proportional to Hp' in Eq. (3.52), the dispersion relation becomes 



where 



Q fu-u' 



Q' 



(3.55) 



(3.56) 



Of course, when the neutrino gas is not moving relative to the matter background ([/' = 0), ^ = 1 and Eq. (3.55) 
reduces to the form given in Ref. . 

The salient feature here is that, in general, the dispersion relation is not isotropic, so that the splitting between 
the two circularly polarized modes is different depending on the direction of propagation of the photon relative to the 
velocity of the neutrino gas. To assess the consequences that this effect can have on the analysis given in Ref. ||l^ we 
consider two extreme cases. 

a. Q perpendicular to U' . Using Q • [/' = 0, it is very simple to verify that 



u ■ u 



while Q' = V[7'2 + Using = -Q'^/q'^, it then follows that 



qVi + if' 
Vn^u'^ + < 



(3.57) 



For small velocities of the neutrino gas this reduces to 1, as it should be, while for large velocities it implies that the 
effect of the Hp ^ term is reduced by the factor Q/fl for fl > Q. 
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b. Q parallel to U' . We set 



Q - XQU' 



(3.58) 



to include the possibility that the photon propagates antiparallel to the velocity of the neutrino gas. A little algebra 
shows that in this case 

, , -u"^Q'' + xnQ\u'\ 



while Q' 



n\u'\-\QU' 
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Therefore 



Q - X(3'fl 

m'-xQ\ ' 



(3.59) 



where we have defined the speed of the neutrino gas 



P' 



U' 

JJIO 



(3.60) 



Eq. ( 3.59 ) reveals in a simple way the anisotropic nature of the dispersion relation. For example, if the velocity of 
the neutrino stream is such that 



(3.61) 



then ^ = — 1 or +1 depending on whether the photon is propagating parallel or antiparallel to /?', respectively. In 
particular, this implies that the frequency difference between the dispersion relations of the two (circularly polarized) 
transverse modes is the opposite to what it is if the neutrino gas is not movi ng relative to the matter background. 
This effect is easy to understand by noticing that, if the condition in Eq. (3.61) holds, then a photon moving parallel 
to (3' appears to be moving in the opposite direction in the rest frame of the neutrino gas. 



IV. OUTLOOK 



Although our work has been largely motivated by the study of the electromagnetic properties of a neutrino gas 
that moves, as a whole, relative to a plasma, the approach is applicable to a wider class of problem in similar physical 
systems, involving relativistic plasmas or high temperature gauge theories. The field theory methods employed here 
allow us to consider situations for which the semiclassical approaches, such as those based on kinetic theory, are 
not suitable, and those for which the full power of the techniques and methods that have been developed for high 
temperature field theory calculations must be employed. Some possible extensions of the present work along those 
lines would involve the calculation of the imaginary part of the self-energy to determine the damping rates, and the 
application of the resummation methods | p3| to study the dispersion relations in those circumstances in which the 
perturbative approximations are not reliable. 
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